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We consider interaction of a femtosecond light pulse with a one-dimensional photonic crystal with 
relaxing cubic nonlinearity in the regime of self-trapping. By use of numerical simulations, it is 
shown that, under certain conditions, the spectra of reflected and transmitted light possess the 
properties of narrow-band (quasi-monochromatic) or wide-band (continuum-like) radiation. It is 
remarkable that these spectral features appear due to a significant frequency shift and occur inside 
a photonic band gap of the structure under investigation. 
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I. INTRODUCTION 

The need for taking into account the noninstantaneous- 
ness of nonlinear response of a medium was realized soon 
after the rise of nonlinear optics. From the end of the 
1960s specialists studied the influence of nonlinearity re- 
laxation in the framework of the Debye model on such 
effects as laser beam self-focusing [l|, 0] and parametric 
amplification Q. Among recent studies, attention has 
been attracted to modulational instability effects in me- 
dia with noninstantaneous nonlinearity [4i-|6|, resulting in 
generation of pulse trains [7,] and solitons 0], instability 
of speckle patterns Q , and reshaping of solitary pulses 

However, the theoretical nonlinear optics of photonic 
band gap materials usually deals with instantaneous pro- 
cesses of nonlinearity. Many results may be found in re- 
views and monographs (see, for example, Refs. pdl - fl^ ). 
We should also note some effects connected with ultra- 
short pulse interaction with nonlinear photonic crystals, 
such as pulse compression and temporal soliton formation 
[Tsl ilGj , subdiffractive propag ation (TtI . [Tsj , and pulse lo- 
calization on a defect |19l - l2]| . 

In this paper we consider spectral transformations of 
femtosecond pulses interacting with a one-dimensional 
photonic crystal with relaxing cubic nonlinearity. As 
was shown in our previous publication f22| . light self- 
trapping occurs in such a nonlinear structure due to for- 
mation of a nonlinear dynamical cavity (or trap) inside 
it. The present paper is a logical continuation of that 
paper. The importance of spectral investigation is con- 
nected with the possibility of spectral broadening, which 
in some extremal cases can result in supercontinuum gen- 
eration. This phenomenon can be observed, for example, 
in photonic crystal fibe rs |23| or in filamentation pro- 
cesses in bulk materials [245. Spectral broadening is one 
of the main points of our research. 

The paper is logically divided into several sections. In 
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Sec. |TT] we give the problem formulation and consider 
some additional details of the self-trapping effect impor- 
tant for the present paper. Section Hill is devoted to the 
spectral features connected with the nonlinear interac- 
tion of a pulse with a photonic crystal in the regime of 
self-trapping. Finally, Sec. IIVI contains a short conclu- 
sion. 



II. SELF-TRAPPING EFFECT 

Propagation of an ultrashort pulse in a one- 
dimensional nonlinear photonic crystal [a structure of 
(A_B)^ type] is described by the Maxwell wave equation 

with the dependence of refractive index on light intensity 
/ = as follows, 

n{z,t) = no{z) + Sn{I,t). (2) 

Here E is the electric field strength, no{z) is the linear 
part of the refractive index varying along the z axis, and 
6n is the nonlinear part of the refractive index, which 
is g overned by the Debye model of relaxing nonlinearity 

dSn 

tnl-^ + on = 7121 , (3) 

where 712 is the Kerr nonlinear coefficient and i„; is the 
relaxation time, which is assumed to be of the order 
of several femtoseconds (fast electronic cubic nonlinear- 
ity). Further we consider femtosecond light pulses with 
the amplitude of Gaussian shape A = exp(— t2/2ip), 
where tp is the pulse duration. To analyze the interac- 
tion of such a pulse with a nonlinear photonic crystal, 
we use the finite-difference time-domain method of nu- 
merical simulations which was described in detail in Rcf. 
[2^ 1. The spectra of pulses (incident, reflected, transmit- 
ted) in this paper are calculated as the absolute values 
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of the Fourier transform of the corresponding field pro- 
files. The spectra are normalized to the peak value of the 
incident pulse spectrum which is recognized as unity. 

The parameters used in our calculations are as follows: 
the linear parts of the refractive indices of the layers A 
and B of the photonic crystal Ua — 2 and rtf, = 1.5, re- 
spectively; their thicknesses a = 0.4 and b = 0.24 /im; the 
number of layers N — 200; the pulse duration tp — 30 
fs; the central wavelength of the initial pulse spectrum 
is Ac = 1.064 /im if not stated otherwise. The nonlinear 
coefficient of the material is defined through the nonlin- 
ear term of the refractive index, so that 712/0 = 0.005; 
this means that the pulse amplitude is normalized by the 
value ^0 — V^o ■ The relaxation time of the nonlinearity 
of both layers is t„; — 10 fs. 

As it was predicted in our previous work p3 [. the in- 
teraction of the pulse (whose duration is comparable to 
the relaxation time) with a nonlinear photonic crystal 
results in the eflfect of pulse self-trapping. This situa- 
tion when the energy of radiation leaving the structure 
is only a small fraction of the incident pulse energy is 
shown in Fig. [TJ It is seen that, for large enough in- 
tensity of the pulse, the output energy demonstrates a 
profound decrease corresponding to the self-trapping of 
the pulse inside the photonic crystal. The output energy 
is calculated by intensity integration over time at the in- 
put and output points of the structure; thus we obtain 
the energies of the reflected and transmitted light, and 
the total output energy is their sum. As the amplitude 
of the incident pulse increases further, light is trapped 
closer and closer to the input face of the crystal so that 
the reflected radiation energy gets larger and larger. Fi- 
nally, the trapping occurs near the very input, so that 
most of the light is immediately reflected. The range of 
pulse durations and relaxation times for the self-trapping 
effect to be observed in our structure was studied in Ref. 
[2^ as well: tni varies from a fraction of a femtosecond 
to about 150 fs, and tp from about 10 fs to about 200 fs. 

For better understanding of this effect, let us consider 
its frequency dependence. In Fig. [2] we see the reflectiv- 
ity and group velocity dispersion (GVD) of the photonic 
crystal considered as functions of the light wavelength. 
It is well known that the GVD parameter k2 = (Pk/dw^ 
is decisive in observation of pulse compression 0, [2^ . 
Indeed, if the nonlinearity coefficient 712 is positive (this 
is the case in our consideration), then one needs to have 
a medium with negative GVD. Since the self-trapping ef- 
fect is characteristic for the regime of pulse compression 
[225 , we can expect that the pulse will be trapped inside 
the photonic crystal if the pulse spectrum lies in the neg- 
ative dispersion domain. This expectation is justifled in 
Fig. [3] where the output energy dependence on the cen- 
tral wavelength of the pulse spectrum is represented. The 
dip in this dependence is unambiguously correlated with 
the negative-GVD region in Fig. [2jb). Moreover, the ap- 
pearance of a minimum in the output energy implies that 
there are some competitive processes which come in con- 
tact and determine the result of the pulse-crystal inter- 
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FIG. 1. (Color online) Dependence of the output light energy 
(normalized to the input energy) on the peak amplitude of 
the incident pulse. Energy is integrated over the time 200tp 
(about six times larger than the pulse transmittance time in 
the linear regime). 
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FIG. 2. Spectral curves for reflectivity (upper panel) and 
group velocity dispersion (lower panel) of the photonic crys- 
tal under consideration. The parameters of the structure are 
given in the text. 



action. It seems natural to suggest that these processes 
are dispersion spreading and nonlinear light-matter in- 
teraction as in the case of usual pulse compression. 

Figure [SJa) also shows that, for larger values of the 
pulse amplitude, the optimal value of the central wave- 
length is situated further from the band gap (compare 
the curves at Am — 2.5Ao and 5Ao). At the same time, 
the dip for Am — 2.5^o is deeper, which is in accor- 
dance with the assumption about optimal (or close to 
optimal) trapping in this case (see Fig. [T]) . Note that for 
Am = 2^0 there is no self-trapping behavior but, nev- 
ertheless, there is a very shallow dip even closer to the 
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FIG. 3. (Color online) Dependence of the output light energy 
(normalized to the input energy) on the central wavelength 
of the incident pulse, (a) The results for different values of 
the amplitude, (b) Comparison of transmitted and reflected 
energies for the pulse with Am = 5Ao. Energy is integrated 
over the time 200t„. 



forbidden gap than at Am = 2.5Aq. 



III. SPECTRAL FEATURES OF LIGHT IN THE 
SELF-TRAPPING REGIME 

The change in shape of the reflected intensity con- 
nected with the dynamical trap formation can be traced 
in Fig. m where we plotted the results for different values 
of incident pulse amplitude corresponding to the cases 
marked by the arrows in Fig. [TJ Most of all, we are 
interested in Figs, ^c) andlH^d), which show some re- 
markable features to be discussed in detail. In contrast 
to the usual peaks of reflected light in Figs. lU^a) and 
Sljb), corresponding to the cases of zero trapping and 
maximal trapping, respectively. Fig. \^c) was calculated 
for a point on the upward slope of the reflected energy 
curve (the amplitude is A„i = 5Aq). This means that 
self-trapping still exists but some part of the radiation 
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FIG. 4. Shape of the reflected intensity at different peak am- 
plitudes of the incident pulse: (a) Am = Aq, (b) Am = 3Aq, 
(c) Am ~ 5Ao, and (d) Am = 7Ao. The central wavelength 
of the pulse spectrum is Ac = 1.064 ^m. 
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FIG. 5. Spectra of (a) reflected and (b) transmitted radiation 
at different peak amplitudes Am of the incident pulse corre- 
sponding to those used in Fig. U The dashed curve depicts 
the spectrum of the incident Gaussian pulse. The band gap 
of the linear photonic crystal is shaded. 
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FIG. 6. (Color online) (a) Shape of the transmitted intensity, 
(b) Spectra of reflected and transmitted radiation. The peak 
amplitude of the incident pulse is Am = 3^o- Calculations 
were performed for the photonic crystal with linear A-layers 
(first layers of the period). The band gap of the linear pho- 
tonic crystal is shaded. 



leaves the nonlinear trap inside the photonic crystal. As 
one can see in Fig. HJ^c), this leaving radiation represents 
almost stationary radiation for quite a long time. As a 
result, we can expect that the spectrum of reflected light 
has to possess a pronounced narrow peak corresponding 
to this quasi-monochromatic radiation. This expectation 
is entirely justified, as spectral plots demonstrate in Fig. 
[Sfa). It is seen that, for Am = 5Ao, the spectrum of 
reflected light really has a sharp peak, while in the cases 
Am ~ Aq and Am = SAq the spectra approximately cor- 
respond to the spectrum of the incident Gaussian pulse. 
Moreover, the position of the spectral peak also provokes 
our interest, since it is situated deep inside the photonic 
band gap of the structure. 

If we take a pulse with greater amplitude {Am — 7Ao), 
then, as mentioned above, the reflected light appears im- 
mediately with a wide and seemingly unstructured enve- 
lope [Fig. m^d)]. It turns out that its spectrum in this 
case completely covers the band gap in continuum-like 
fashion [see Fig. [5ja)]. Note that the spectra of trans- 
mitted radiation do not intrude into the forbidden gap, as 
is witnessed by Fig. ^h). In general, spectral broaden- 
ing can be linked with self-phase-modulation resulting in 
the generation of new frequencies in the pulse spectrum 
due to the temporal variation of the refractive index [^J] . 
However, the fact that the sharp edges of the spectrum 
include just the entire band gap seems to be unexpected. 
In fact, we have a situation when light converts under 
nonlinear interaction in such a way that the spectrum is 
more and more pulled into the forbidden gap. On the 
other hand, we should keep in mind that self-trapping is 
connected with a local change of reflective properties of 
the photonic crystal 22|. 

Thus, if we are on the upward slope of Fig. [1] we 
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FIG. 7. (Color online) Spectra of reflected and transmitted 
radiation at different central wavelengths Ac. The peak ampli- 
tude of the incident pulse is Am ~ 5Ao. The band gap of the 
linear photonic crystal is shaded. The position of the spec- 
trum of the incident pulse is characterized by the bell-shaped 
curve in the spectrum of the reflected light. 



can obtain narrow-band radiation in reflection. In other 
words, this corresponds to stronger coupling between the 
pulse and the nonlinear structure than in the case of the 
optimal self-trapping effect due to the greater value of 
the incident intensity. One might suggest that some- 
thing similar should be observed in the opposite situa- 
tion when the light-medium interaction gets weaker, i.e. 
in the region of the very abrupt downward slope in Fig. 
[2 Obviously, the narrow-band spectrum is expected to 
be obtained in transmitted (not reflected) radiation in 
this case. In order to prove this statement, we use an- 
other method to make the coupling between the pulse 
and the nonlinear photonic crystal weaker. We take ma- 
terials with smaller nonlinearity, rather than decreasing 
the intensity of the pulse. In Fig. [6]the results are shown 
for the (yli3)^-structure with linear A layers, while the 
parameters of the B layers remain unchanged. It is seen 
that the structure of the transmitted radiation is simi- 
lar to that of Fig. [D^c). As a result, in the spectrum 
of transmitted light a pronounced quasi-monochromatic 
peak occurs [Fig. Elb)]. This peak, however, is situated 
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FIG. 8. (a)-(c) Spectrograms of reflected radiation at peak amplitudes of the incident pulse Am = ^o, 5Aq, and 7Ao, respectively. 
The central wavelength is Ac = 1.064 fj,m. (d) Spectrogram of transmitted radiation at Am = 5Ao and Ac = 0.8 fj,m. 



near the very edge of the band gap of the hnear photonic 
crystal. Obviously, low-intensity quasi-monochromatic 
radiation seen in Fig. [BJa) cannot significantly change 
the refractive properties of the structure through which 
it is to be transmitted. Therefore, there is only a slight 
shift of the forbidden gap, which can be referred to as a 
self-induced transparency effect in the nonlinear photonic 
crystal. 

Let us return to Fig. [3] and consider spectral transfor- 
mations for incident pulses with different central wave- 
lengths Ac. The dip corresponding to the self-trapping 
phenomenon is situated between two band gaps plotted 
in Fig. [5] (we call it the inter-gap region). Changing Ac 
in this region, one can obtain all the variants of spectral 
peculiarities discussed above and even more as can be 
seen in Fig. [71 Further, we list the main features seen in 
this figure; 

(i) When the spectrum of the incident pulse is out of 
the inter-gap region (Ac — 0.71 /im, which is the region 
of negative GVD; see Fig. [2]), the spectrum of reflected 



radiation has a sharp peak near the very low-frequency 
(on the dip side) edge of a narrow band gap. 

(ii) If the spectrum is in the positive-GVD domain 
(Ac = 0.8 /Ltm), the peak for reflected radiation still oc- 
curs, but there also appears a peak in the transmitted 
light spectrum near the edge of the wider (long-wave) 
band gap. Self-induced transparency is also observed 
due to the shift of the forbidden gap. Obviously, this 
situation corresponds to the weak light-matter coupling 
regime discussed previously in connection with Fig. [51 

(iii) As we move further inside the negative-GVD re- 
gion (the region of the self-trapping dip, Ac = 0.9 and 
1 /.tm), the peaks near both the wide and narrow gaps 
diminish and become less and less pronounced (the peak 
for transmitted light moves away from the edge of the 
gap). 

(iv) At Ac = 1.06 ^m we see a narrow quasi- 
monochromatic peak actually inside the forbidden gap 
(compare with Fig. [51 Am = 5Aq). Note that the re- 
flected radiation also appears inside the narrow (short- 
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wave) band gap. 

(v) At Ac = 1.08 /iin the reflected light spectrum in 
the band gap widens and takes a continuum-like shape 
with characteristic oscillatory fine structure (compare 
with Fig. [S]at A™ = 7Ao). 

(vi) If the initial spectrum is almost entirely inside the 
band gap (Ac = 1.15 /im), then we have the usual reflec- 
tion without any evidence of nonlinear interaction. 

(vii) At Ac = 1.4 /im we are also outside the dip and 
inside the positive-GVD region. The spectra of both re- 
flected and transmitted light demonstrate a sharp break 
right on the edge (the side opposite to the dip) of the 
band gap. There is no any sign of self-induced trans- 
parency. 

Finally, to make clear the connection between tem- 
poral curves and spectra, we turn to the spectrogram 
technique widely used in supercontinua investigation [23[ . 
The spectrogram is calculated as [1^ 



+ 00 



E{t)g{t- 



(4) 



where E{t) is the field under investigation (in our case, 
the reflected or transmitted field), g{t) = exp(— t^/2tp) 
is the gate function which is chosen to be a replica of 
the input pulse. The spectrogram S{uj, r) allows an intu- 
itive understanding of the correlation between temporal 
and spectral features of a given signal. In Fig. [S] such 
spectrograms are shown; Figs. [8Ua)-[8l[c) corresponding 
to temporal and spectral curves of the reflected radia- 
tion depicted in Figs. |l]and[5la), respectively. Figure 
ISKd) represents the spectrogram of transmitted radiation 
of the spectrum demonstrated in the upper right panel 
of Fig. [7] (the case of Ac = 0.8 ^m). 

The spectrogram of reflected light at the incident pulse 
amplitude A„i — Aq [Fig. UKa)] shows two intensity peaks 
seen in Fig. IDJa). These peaks are concentrated near 
the central wavelength Ac — 1.064 /xm without any sig- 
nificant frequency shift. Such a shift is easily seen in 
Fig. [S^b) at Am = 5Ao, so that low-intensity quasi- 
monochromatic radiation occurs exactly inside the pho- 
tonic band gap. It is also worth noting that at t > 75tp 
a frequency shift in the reverse direction (toward Ac) ex- 
ists. This fact can be associated with the chaotic ending 
of quasi- monochromatic radiation seen in Fig. Hfc). For 
larger input intensity (A,„ = 7Ao), this reverse shift oc- 
curs earlier in time, but the spectrum covers a wider fre- 
quency range with approximately uniform intensity giv- 
ing rise to continuum-like radiation inside the band gap. 
The last spectrogram [Fig. [S^d)] shows the transmitted 
quasi-monochromatic radiation with a large frequency 
shift from Ac = 0.8 fxm to about 1.12 /im. It is easily 
seen that this shift happens very fast in time. 



IV. CONCLUSION 

In conclusion, in this paper we have studied the spec- 
tral transformations of ultrashort (femtosecond) light 
pulses resulting from their interaction with a nonlinear 
photonic crystal in the regime of self-trapping. In our 
analysis we used only the processes of light self-action, 
so that the effect of generation of wide and narrow spec- 
tra cannot be connected by high harmonics and sum- 
frequency appearance. However, as our results demon- 
strate, these self-interaction processes are sufficient for 
impressive spectral transformations in nonlinear photonic 
crystals. These transformations concern both reflected 
and transmitted light spectra and depend on the regime 
of light-material interaction. In particular, if this inter- 
action is strong (the pulse is trapped near the entrance of 
the photonic crystal) , a narrow peak and continuum-like 
spectral features occur in reflected light. On the other 
hand, if the light-structure interaction is weak (the pulse 
is trapped near the exit of the photonic crystal), a nar- 
row peak near the edge of the band gap appears in the 
transmitted light spectrum. Obviously, relaxing (nonin- 
stantaneous) behavior of the nonlinearity and periodic 
change of the linear refractive index (a photonic crystal 
per se) are the key conditions due to the necessity of 
self-trapping. 

We should also say a few words about perspectives of 
this research. First, some improvements are possible in 
the realization and control of the self-trapping effect by 
adjustment of the photonic structure. In particular, a 
chirped photonic crystal with varying period can be em- 
ployed to shift the trapping position inside the structure 
and, perhaps, to relax the requirements on the materi- 
als. However, this modification is still to be studied in 
detail. The second question is connected with the possi- 
bility of experimental realization of self-trapping and the 
corresponding spectral effects. Although the parameters 
used do not belong to some specific nonlinear medium, 
they seem to be quite realistic. The relaxation times 
(a few femtoseconds) are characteristic of media with a 
fast electronic mechanism of Kerr nonlinearity. How- 
ever, such media possess relatively low nonlinear coef- 
ficients; therefore one has to use high-intensity pulses 
(~ 100 GW/cm^) and take into account the problem of 
the damage threshold, which is high enough in the case 
of femtosecond pulses. We believe that one can find some 
materials (for example, doped glasses [2^]) which satisfy 
all these conditions. 
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